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We investigate the effects of a temperature-dependent shear viscosity 
over entropy density ratio 77/s on the flow anisotropy coefficients V2 and Vi 
in ultrarelativistic heavy- ion collisions at RHIC and LHC. We find that V4 is 
more sensitive to the viscosity at low temperatures than V2 ■ At RHIC V2 is 
mostly affected by the viscosity around the phase transition, but the larger 
the collision energy, the more the quark-gluon plasma viscosity affects V2 ■ 

PACS numbers: 25.75.Ld, 12.38.Mh, 24.10.Nz 

1. Introduction 

Presently, most works aiming at the determination of the shear viscosity 
of strongly interacting matter assume a constant shear viscosity over entropy 
density ratio, 7]/s. However, this ratio can be a strongly varying function 
of temperature both in hadronic matter and in the quark-gluon plasma. In 
this work we study consequences of such a temperature dependence [1]. 

We model the space-time evolution of matter formed in heavy-ion col- 
lisions using relativistic dissipative hydrodynamics [2]. We assume longitu- 
dinal boost invariance and neglect the net-baryon number. Essential inputs 
to the model are the equation of state, the initial state and the transport 
coefficients. We consider here only the shear viscosity. 

As equation of state we use a recent lattice parametrisation [3] with 
chemical freeze-out at T = 150 MeV. The initial energy density at tq = 
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Fig. 1. rj/s parametrizations. 



1.0 fm is proportional to the density of binary nucleon-nucleon collisions in 
the transverse plane. The maximum energy density is fixed to reproduce the 
measured multiplicity in the most central collisions [HIS]. For s/s^n = 5.5 
TeV Pb+Pb collisions we use the multiplicity predicted by the minijet + 
saturation model [B]. To compensate for different entropy production for 
different parametrizations of the shear viscosity, the initial energy-density 
profiles are normalised differently for each parametrization. Freeze-out is 
implemented using the Cooper-Frye formula [7] on a T^ec = 100 MeV hyper- 
surface including the dissipative correction 5f to the thermal distributions. 

For r]/s, we consider the four different parametrizations shown in Fig. [TJ 
The minimum value of ij/s is fixed to he r]/s = 0.08 at T = 180 MeV for 
all parametrizations. The shear relaxation time [2j is taken to be tr = 
5T]/{e+ p), where e and p are energy density and pressure, respectively. 

2. Results 

The elliptic flow coefficient, V2{pt), for charged hadrons in the 20-30% 
most central collisions at RHIC, at the present LHC energy, and at the full 
LHC collision energy is shown in the left panels of Figs. [2HH respectively. 
We note that at RHIC, the high-temperature part of r]/s has practically no 
effect on the results. On the other hand, the viscous suppression of V2{pt) 
is strongly enhanced if we increase the hadronic r//s. In low-energy LHC 
collisions, both hadronic and QGP viscosity have a similar effect, whereas 
at the full LHC energy the behaviour is opposite to that seen at RHIC: 
V2{pt) is almost independent of the hadronic r//s, but sensitive to the high- 
temperature viscosity. 

The anisotropy coefficient V4^{pt) (right panels of Figs. [2H1]) exhibits a 
similar pattern where the sensitivity to the low-temperature viscosity de- 
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Fig. 2. V2(pt) (left) and V4{pt) (right) of charged hadrons in the 20-30% most 
central Au+Au coUisions at ^/snn = 200 GeV (RHIC). Data are from Refs. 



creases, and the sensitivity to the high-temperature viscosity increases with 
increasing collision energy. However, in general U4(pr) is more sensitive to 
the shear viscosity in the low-temperature region than V2{pt)- At the low 
LHC energy V4{pt) behaves like V2{pt) at RHIC, and at the full LHC energy 
V4:{pt) behaves like V2{pt) at low LHC energy. 

To further investigate the sensitivity of V2 and f 4 to viscosity at different 
temperatures, we devised a set of simple parametrizations for r]/s{T), see 
Fig. [5l We take r]/s = 0.08, except in the vicinity of a temperature Tp^ak, 
where ?7/s(rpeak) = 0.24, and the width of the peak is 10 MeV. By varying 
^peaki we can study at what temperature r]/s affects V2 and 1^4 most. 

The results are shown in Fig. [6] as the relative change with respect to u 2 
and Vi evaluated using a constant rj/s = 0.08 during the entire evolution. 
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Fig. 3. V2{pt) (left) and V4{pt) (right) of charged hadrons in the 20-30% most 
central Pb-^Pb collisions at ^/sWn = 2.76 TeV (LHC). Data are from Ref. [TUl. 
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Fig. 4. V2{pt) (left) and V4{pt) (right) of charged hadrons in the 20-30% most 
central Pb+Pb coUisions at ^/sWn = 5.5 TeV (LHC). 



It is seen that at RHIC, the viscosity around T = 170 MeV has the largest 
effect on V2, but the region of strongest sensitivity moves to larger temper- 
atures and becomes wider with increasing collision energy. If we ignore the 
point at T = 110 MeV (because its main effect is via 6f), the behaviour 
of V4 is slightly different: With increasing collision energy the suppression 
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at large temperatures increases, and at low temperatures decreases, but the 
temperature where r]/ s suppresses most hardly changes. 

An interesting feature is that there is a region where r]/ s does not sup- 
press flow anisotropics, but the larger the r]/s the larger the anisotropy! 
This can be understood in the following way: At early times the main ef- 
fect of shear viscosity is to inhibit the longitudinal expansion, and enhance 
the transverse expansion, instead of reducing the difference between the ex- 
pansion in in-plane and out-of-plane directions. Thus, shear viscosity leads 
to larger transverse flow velocity. A simple blast wave model [TT] demon- 
strates that if nothing else changes, a larger transverse flow velocity leads 
to a larger V2{pt) of light particles. At RHIC a similar reasoning leads to 
the insensitivity to the plasma viscosity, since the effects of increasing flow 
velocity and smaller difference between in-plane and out-of-plane directions 
cancel each other. 

All this does not mean that V2 would not be formed early. To take 
into account the thermal motion during the evolution, we characterise the 
time-evolution of V2 by evaluating the V2 of fictitious m = 140 MeV bosons 
at different times Tj. We use the Cooper-Frye formula on a hypersurface 
consisting of two parts: A constant temperature hypersurface with T = 
Tdec = 100 MeV and t < Ti, and a constant time hypersurface with t = ti 
and T > 100 MeV. This approach has the advantage that at the end of the 
evolution it matches the V2 of thermal pions without any adjustment. 

As seen in Fig. [71 V2 is built up early, but the effect of rj/s at early 
times is relatively small. The left panel depicts the evolution at RHIC, and 
one can see that at r = 2 fm, the larger rj/s above the transition region 
has reduced the V2 somewhat, but this difference is soon erased and the 
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Fig. 7. Time-evolution of the pT-averaged V2 of fictitious m = 140 MeV bosons 
("tt") in 20-30% most central Au-|-Au (left) and Pb-|-Pb (right) collisions at 
^5]^ = 200 GeV (RHIC) and .JsWn = 2.76 TeV (LHC), respectively. 
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low-temperature viscosity dominates the evolution from r ~ 3 f m on. The 
situation is different in the evolution at the lower LHC energy depicted 
in the right panel. Now the plasma viscosity causes a clear difference at 
r ~ 2 fm, and this difference persists to much later times than at RHIC. 
Eventually the hadronic viscosity takes over and reorders the curves, but it 
cannot completely compensate the differences built up earlier. 

In summary, we have investigated how the temperature dependence of 
rj/s affects the flow anisotropy coefficients V2 and v/^. We found that the 
temperature where Vn is most sensitive to viscosity varies with the collision 
energy — the larger the collision energy, the larger the temperature where 
the suppression is strongest. We also saw that v^^ is most sensitive to the 
viscosity at a lower temperature than V2- It remains to be seen whether this 
is a general trend: the larger the n, the lower the temperature where Vn is 
most sensitive to the value of r//s. 
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